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Abstract: This paper deals with the simulation of a permanent magnet synchronous machine (PMSM) 
to be developed for driving a scroll compressor of a vehicle air conditioning system. The simulations 
were needed for the verification of the pre-sizing and for the electromagnetic behavior analysis 
of the PMSM. These were performed by using the ANSYS Motor-CAD software. All the obtained 
results emphasized that the pre-sizing of the PMSM was performed correctly, and the designed 
electrical machine can fulfill all the prescribed requirements. To highlight the superiority of the 
chosen surface-mounted PMSM topology, a comparative study was performed for four PMSMs 
having the same stator but different rotor topologies having the same quantity of permanent magnets. 
The study revealed that the PMSM with permanent magnets on its rotor surface is the best-fitting 
variant for the foreseen application. By intensively using the ANSYS Motor-CAD software, the authors 
appreciated that this product is a very useful tool in the hands of the designers not only in evaluating 
the machine’s computed geometry (by checking its correctness and feasibility), but also in its in-depth 
electromagnetic analysis and in calculating its parameters of interest. 


Keywords: electrical machine design; permanent magnet synchronous machine; design software; 
ANSYS Motor-CAD program; finite element method; numerical electromagnetic field computation 


1. Introduction 


Even though electrical machines are regarded as mature products, their design is 
a complex, time-consuming, and demanding engineering task as it requires profound 
multidisciplinary knowledge of electrical engineering (electrical machines and drives, elec- 
tromagnetic field theory, electrical materials), mathematics (numeric methods), mechanical 
engineering, fluid mechanics and thermodynamics, acoustics, manufacturing technology, 
technical drawing, computer sciences, etc. [1]. Additionally, the electrical machine design- 
ers must be knowledgeable not only about the given specifications (which refers to machine 
ratings, customer performance expectations, etc.), but also of the existing standards, labor 
and materials costs, manufacturing requirements and restrictions, etc. 

All electrical machine designers are facing a great variety of issues, some of which may 
not have a single answer, but rather several. The computations needed to be performed 
during the sizing of the electrical machines are complex and almost repetitive due to the 
frequent recalculations because of the non-fulfillment of the imposed conditions. Therefore, 
it is mandatory to apply some computer programs to perform the needed computations. 
Moreover, nowadays, the sizing of an electrical machine must be completed by adequate 
optimization to achieve the best structure within certain constraints [2,3]. 

Nowadays, it is not sufficient that the developed electrical machine works well. Global 
attention has been focused on the rising cost of energy over the past years, and it appears 
that this trend will continue. As a result, it is now more crucial than ever to design 
electrical machines with very high efficiency. Their use can significantly reduce world 
energy consumption [4]. In the post-pandemic period, during wartime, all the industries 
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are facing massive unprecedented supply chain disruptions [5]. While the supply chain 
management strategies are re-examined, the material savings and recycling possibilities 
must also be the focus of the designers of any products [6]. 

Electrical machine designers have at their disposal numerous software tools to ease and 
improve their challenging engineering work. Most of them are using professional finite ele- 
ment analysis (FEA) software products for numeric electromagnetic field computations [7,8], 
such as JMAG-Designer [9,10], ANSYS Maxwell [11,12], Simcenter Magnet [13,14], COMSOL 
Multiphysics [15,16], Flux 2D [17-19] or its three-dimensional variant Flux 3D [20,21]. All 
such programs can perform precise computations for obtaining the magnetic field distri- 
bution within the electric machine, the generated forces, magnetic, material-related, and 
dynamical parameters, etc. They have versatile graphical user interfaces (GUIs) helping the 
engineers in adapting their models to the specific needs of the design. The keys to obtaining 
very precise results are the applied advanced adaptive mesh generators [22,23] and the 
solvers based on the most effective, time-saving, and accurate numerical methods [24,25]. 
Their embedded scripting tools enable the designer to simulate a great diversity of electrical 
machine topologies by only changing some basic parameters of the script and thus totally 
controlling the preprocessing phase of the simulations. 

All commercial FEA software is interoperable, enabling the use of different external 
tools offered by other programs. This is performed by flexible communication with the 
other software products. The link of the FEA programs to Simulink must be primarily 
mentioned, which permits the easy development of the electrical machines’ control system 
and thus their precise dynamic analysis [26]. These programs can export and import 
graphical files of a great variety of types, too. 

There are specialized software solutions on the market that not only carry out a 
general FEA, but also greatly aid electrical machine designers by providing unique features 
essential in sizing and analyzing such devices. For example, Simcenter SPEED provides 
a variety of machine templates and libraries to assist users in defining the geometry, 
winding, drive, and control strategy of an electrical machine under development to gain 
information from simulations about its current, torque, speed, losses, etc. [27]. Another 
similar software product, Simcenter Motorsolve, has an easy-to-use template-based GUI 
for defining a great diversity of electrical machine topologies, which can be analyzed by 
applying both equivalent magnetic circuit-based calculations and FEA [14,28]. 

A special place in the market of electrical machine design and analysis software is 
occupied by ANSYS Motor-CAD. This useful program can perform in-depth multi-physical 
analysis of a great diversity of electrical machines [1,29,30]. The users have at their disposal 
four distinctive modules [31]. The main component is the electromagnetic one, which is 
designed for the electromagnetic performance assessment of the electrical machines by 
using precise analytical methods and FEA [7,8]. By applying the thermal unit, a 3D lumped 
thermal circuit can be built up automatically, enabling the user to determine the steady 
state and transient thermal characteristics of the machine taken into the study [32-37]. The 
third module, named Lab, is a special toolbox for combined electromagnetic and thermal 
simulations, which can ease the modeling and optimization tasks of the designers [38-41]. 
The mechanical module of the ANSYS Motor-CAD software uses the effective built-in 
FEA solver to calculate miscellaneous mechanical properties. This component can be 
very useful in the electromagnetic Noise, Vibration, and Harshness (NVH) analysis to 
be performed [42,43]. 

Our studies were focused on the development of an electrical machine for the com- 
pressor of the heating, ventilation, and air conditioning (HVAC) system of a car. Nowadays, 
especially in the summer, air conditioning is a vital necessity because it is the only way to 
circulate and cool the air inside any vehicle. Furthermore, these devices also clean and dry 
the indoor air [44,45]. Considering the main advantages of permanent magnet synchronous 
machines (PMSMs), like high efficiency, great specific power and torque, and good behavior 
at high speeds, such a machine was proposed for the given specific automotive application. 
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The paper begins with a brief survey of the software tools used by electrical ma- 
chine designers. The focus is set on the ANSYS Motor-CAD software, since in the paper, 
the authors intend to share their experience in using this product in designing PMSMs. 
Next, the scope of designing the certain PMSM and the listing of its established design 
specifications follow. In this part of the paper, the results of the common pre-sizing of 
the machine are also detailed. The most consistent part of the paper is dealing with the 
surface-mounted PMSM model construction in ANSYS Motor-CAD and with detailing the 
obtained simulation results. By using the gained experience in using this software product, 
a comparative study on diverse PMSM topologies was performed. Four PMSMs having the 
same stator and permanent magnet quantities but different rotor structures were compared 
from diverse points of view [46-48]. The main aim of this study was to find out whether the 
other three PMSMs exceed the performances of the initial variant having surface-mounted 
permanent magnets. The study evidenced both the superiority of the surface PMSM variant 
and the practicality of applying ANSYS Motor-CAD software in designing and analyzing 
electrical machines. 


2. The PMSM to Be Studied 
2.1. Automotive Applications and the Electrical Machine Type Selection 

A scroll compressor, similar to that given in Figure 1, was proposed for a car HVAC 
system [49,50]. They have the advantage of containing only two essential pieces (a fixed and 
an orbiting scroll) and thus can operate more quietly and smoothly than the reciprocating 
compressors [51]. Moreover, they are more reliable and energy-efficient and less prone to 
mechanical failure due to their smaller number of parts. 


; a’ 
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Figure 1. The scroll-type compressor. (a) The two co-wound spiral-shaped scrolls; (b) internal view 
of the compressor. Courtesy of the Vacuum Science World [52]. 


After a careful analysis, the PMSM was selected to propel the scroll compressor [53]. 
The main reason was that both efficiency and power density is quite high for these electrical 
machines. Additionally, they require little maintenance effort and produce little heat. Due 
to all of these factors, they are among the most significant industrial competitors of the 
very extensively used induction machines [54]. 

PMSMs are widely used in the automotive industry [55]. In 2021, 90% of the newly 
sold full and hybrid vehicles were powered by PMSMs [56]. Other significant automo- 
tive application areas for PMSMs include electric vehicle auxiliaries, like power brakes 
and steering. 

Their global market share was USD 13.45 billion in 2021, and by 2028 it is anticipated to 
surpass USD 22.6 billion. This should result in a notable (7.57%) increase in the compound 
annual growth rate (CAGR) over this time span [57]. Given that the cost of the rare earth 
metals required to produce permanent magnets has more than doubled year over year, this 
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gain in PMSM market share is particularly impressive and can be explained only by the 
outstanding properties of the PMSMs [56]. 

All of this emphasizes the significance of research to enhance PMSMs through the use 
of cutting-edge design tools, new materials, and novel manufacturing procedures [58]. 


2.2. The Design Specifications of the PMSM to Be Developed 


The PMSMs to be studied here and designed to drive an automotive scroll compressor 
must have the mechanical characteristic provided in Figure 2. 


5 


Torque [N-m] 


0 2000 4000 6000 8000 10000 
Speed [r/min] 


Figure 2. The mechanical characteristic of the PMSM to be designed. 


The main input design data for the PMSM were established upon the torque and speed 
development of the given scroll compressor. During the design procedure, the following 
starting data were used: 


number of phases: 3; 
number of pole pairs: 3; 
DC bus voltage: 72 V; 
rated power: 4 kW; 
rated speed: 8000 r/min; 
rated torque: 4.77 N-m. 


2.3. Pre-sizing the PASM 


The first step in designing the PMSM was the pre-sizing. This consisted of computing 
the initial geometry and winding of the PMSM by using simple analytical equations and 
some general rules [59-61]. During this design stage, the iron core and permanent magnet 
volumes, the number of slots, the main sizes of the stator and rotor, the winding diagram 
and the number of turns, the materials to be used, etc., were estimated by the well-known 
relationships among the main machine sizes and other parameters [62-65]. 

The main dimensions of the pre-sized three-phase, 18-slots surface PMSM are given in 
Figure 3. 
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Figure 3. The main sizes (in mm) of the stator and rotor iron core laminations. 


The cross-section of the rotor with the permanent magnets with its main sizes can be 
seen in Figure 4. 


1.8 


Figure 4. Cross-section of the rotor with the permanent magnets on it (all sizes are in mm). 


In the 18 stator slots, a three-phase, double-layer, distributed winding was proposed 
to be placed. As can be seen in Figure 5, the pitch of the winding is three slots. 


Figure 5. The winding diagram of the PMSM. 


Upon these initial data, the model of the PMSM was built up in ANSYS Motor-CAD 
for evaluating the computed geometry and the selected winding, and for checking their 
correctness and feasibility. 
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3. The ANSYS Motor-CAD Model of the PMSM 


By using the ANSYS Motor-CAD program, the building up of the surface-mounted PMSM 
model was very easy, as it needed only a few simple well-assisted steps to be performed. 

To investigate the electromagnetic behavior of this PMSM, the electromagnetic mod- 
ule (named E-Magnetic) of the software was used. This enables the calculation of the 
electromagnetic performances by using both precise analytical and FEA methods [7,8]. 

The program offers to the designers eight predefined electromagnetic templates to be 
used for the study of the so-called Brushless Permanent-Magnet Motors (BPMs) category. 

The first step in the model implementation was the definition of the stator and rotor 
iron cores. As can be seen in Figure 6, the sizes of the transversal and axial cross-sections of 
the PMSM had to be set from the Geometry menu of the program. 


O Geometry |[] winding | WA input Data | {It} Calculation | GS E-Magnetics | EE] Output Data | [Graphs | $ Forces | BZ Sensitivity | 
O Radial |F aia |020 | 


Slot Type: [Tapered Sot v | Rotor Type: | Surface Breadoaf v 
Stator Ducts: [None ~ | Rotor Ducts: | None 
Stator Parameters | Value || Rotor Parameters | Value 

eno EA- E 

Stator Lam Dia 118 || Magnet Thickness 35 

| Stator Bore | 50 ||Magnet Reduction | 0.3918 
[sot with Bottom) | 10.0809 [Magnet Arc [ED] 1135 
| Slot Width (Top) 5,9434 | Magnet Segments 1 
| Slot Depth 24 || Airgap 0.65 
[Sot Comer Radus | 1 |[Banding Thickness 0 
| Tooth Tip Depth 1 || Shaft Dia 25 
| Slot Opening 2 || Shaft Hole Diameter 0 
| Tooth Tip Angle 0 
| eeve Thickness 0 


Geometry Parameterisation: 
@ Dimensions O Ratios Redraw 


(a) 
OGceometry |[] Winding | WA input Data | Hl} Calculation | Æ E-Magnetics | ZE] Output Data | K Graphs |<) Forces | BZ Sensitivity | 
@ Radai Baca |0020 | 


Shaft Type: | Solid 


Radial Ducts: | None v 


ee N E 
| Stator Lam Dia 118 || Motor Length | “120 
Stator Bore 50 || Stator Lam Length 0 
| Airgap 0,65 _ || Magnet Length 0 
Banding Thickness 0 || Magnet Segments 1 
Sieve Tickness | o _|[Rotor Lam Length 20 
Magnet Thickness 3.5 || EWdg Overhang [F] 30 
| Shaft Dia 25 || EWdg Overhang [R] 30 
Shaft Dia [F] 25 || Wdg Extension [F] 5 
Shaft Dia [R] 25 || Wdg Extension [R] 5 
Shaft Hole Diameter 0 || Shaft Extension [F] 30 
| Shaft Extension [R] 0 


Redraw Draw Cooling KZ Draw datum 


Figure 6. The panels of the ANSYS Motor-CAD program for setting the geometrical data of the stator 
and rotor iron cores. (a) Radial view; (b) axial view. 


In the text boxes of the two panels, the main sizes of the iron cores (exterior and 
interior stator diameter, airgap length, shaft diameter, active length of the machine, and 
details about the stator teeth and slots) must be established. 
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Taking advantage of the possibility to export the cross-section of the iron cores in dxf- 
type CAD files, the resulting and the initial (obtained upon the pre-sizing of the machine) 
drawings were overlayed (see Figure 7). As can be seen, the designed and modeled 
structures are very close, so it can be stated that the iron core structure of the surface PMSM 
in the study was precisely modeled in ANSYS Motor-CAD. 


Figure 7. The overlayed drawings of the designed (black lines) and modelled (lines drawn in red) 
iron cores of the PMSM (showing only half of the motor cross-section). 


The winding diagram of the PMSM given in Figure 5 can also be very easily imple- 
mented in ANSYS Motor-CAD from its Winding menu. The lap type of the winding with a 
central path had to be selected. By specifying only five winding parameters (the numbers 
of the phases, turns, parallel paths, layers, and the throw length in slot steps), its model 
could be simply completed. 

In the Pattern panel of the Winding menu, two types of winding diagrams automatically 
generated by the program can be visualized (see Figures 8 and 9). By using these, the 
modeled winding diagram can be compared easily with that previously set up by the 
machine designer. 


© Geometry []Winding |Q input Data | Hf} Calculation |G E-Magnetics | EE] Output Data | 4 Graphs | EG Senstivty | 
CPattem |E] Definition | 


Phase Distribution: 
Independent aa Tas: 


Nx 3phase Winding Layers: 


O 
f 
g 
7 
3 


@ Radial Patem |f Linear Pattem | Phasors| MMF | Hamonics| Factors} 
AlPhases | Ph1 | Ph2 | Ph3 | 
Slot ae | a 


o o uona w N = 


keas daisi S 
| 

bomookoaonocen 
relelalele|aleyelalelo 
þh m|olo|alololalolo|alo 


Hala 
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Figure 8. The radial pattern of the PMSM winding diagrams implemented in ANSYS Motor-CAD. 
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Figure 9. The linear pattern of the PMSM winding diagram as provided by the ANSYS Motor-CAD 
software. 


In the same panel, the phasor diagram, the MMF distribution in the airgap, and its 
harmonic content can be also visualized and checked. 

The Definition panel of the same Winding menu offers a specific and very useful feature 
of the ANSYS Motor-CAD software product. As can be seen in Figure 10, the filling of the slot 
is visualized upon the imposed number of turns, sizes of the conductors and insulators, etc. 
Additionally, a great variety of winding parameters are shown in this window, such as 
different areas, fill factors, etc. 


otor- 


Output Parameter Value | | _ Output Parameter | Value 
Conductors/Slot Drawn [312 | | Slot Area 185,7_| [Covered Wire Area 89.06 | 
Wire Slot Fill (Wag Area) [0.5516 | [Winding Area (+ Line) [171.6 _| [Copper Area 70,25 | 
Wire Slot Fill (Slot Area) [0.4795 | [Winding Area 161.5 _| |impreg Area za | 
Copper Slot Fill (Slot Area) |0.3782 | | Winding Depth 20,518 | [Wedge Area 14.18 
Heavy Build Slot Fill 0,7023 


Figure 10. The Winding-Definition panel showing the way the conductors of the windings are placed 
in the stator slots and listing numerous parameters of the winding. 


The Conductors/Slot Drawn output parameter value should be colored in red to indicate 
an error if the designed winding does not fit into the slot area with the specified winding 
settings at the required fill factor. 

In the next step of the modeling, the materials used in the construction of the surface 
PMSM had to be established. For this purpose, ANSYS Motor-CAD offers a huge database 
of the usual materials of electrical machines. 

Fortunately, all the materials selected during the pre-sizing of the PMSM could be 
found in the provided database: M470-50A-type non-oriented 0.5 mm thick electric steel for 
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the iron cores [66], pure copper for the windings, and the N35UH-type NdFeB permanent 
magnet [67], as can be seen in the material selection panel provided in Figure 11. In ANSYS 
Motor-CAD, there exists the possibility for the users to create their own databases for those 
materials which had not been included in the program database. 


© Geomety | F] winding input Data | {1} calculation | @% E-Magnetics | EE] Output Data | K Graphs | 5% Sensitivity | 
Bia Materals | Settings | & Material database | 
Component | Material from Database | Hectrical | Temp Coef | Magnet | Magnet | Temp | Density | Weight Notes a 
Resistivity | Hectrical | Brat | Relative | Coef Br 
Resistivity | 20°C | Permeablility 
| Units Ohm. | Tesla aC | kg/m kg 
| Stator Lam (Back Iron) M470-50A vl 3.6607 o | | | 7650 | 2.266 
E Lam (Tooth) M470-50A 3.6E07 0 | | 7650 147 
| Stator Lamination [Total] | | | | 3,736 
| Amature Winding [Active] | Copper (Pure) ~| 1.724E08 | 0,003862 | | 8933 1.017 
| Amature EWdg [Front] Copper (Pure) 1.724E08 | 0.003862 | | 8933 | 0.4147 
| Amature EWdg [Rear] Copper (Pure) 1,724E-08 | 0.003862 | | 8933 | 0.4147 
[Amature Winding [Total] __ aes |e | | | 1846 _ 
| Slot Wedge v o 0 1000 | 0.02298 
| Rotor Lam (Back Iron) M470-50A v| 3.6607 0 | 7650 | 0.6537 
| Rotor Lamination [Total] EE | | _ | 0.6537 
[Magnet N35UH v| 1.8606 E 105 | 12 | 70 | 0.1881 
[Shaft [Active] v 0 o. | 7800 | 0.3446 
| Shaft [Front] 3] E7 o po) 7200 | 0.1723 
| Shaft [Rear] v 0 o | 7800 | 0.05743 
[Shaft [Total] | | | | 0.5743 
[Total | | | 7.021 _|Weight [Total] 


Figure 11. The Materials panel of the Input Data menu used for selecting the materials for the different 
components of the modeled PMSM. 


As can be seen in Figure 11, from this panel, information on the main characteristics 
of the used materials can be found. Important data regarding the masses of the PMSM 
components and the total machine mass can also be obtained from here. 

With this last step, the model development in ANSYS Motor-CAD was practically 
finished. Next, the simulation conditions must be set in the Calculation menu of the 
program shown in Figure 12. 


© Geometry | F] Winding | Qf input Data fff} Catcuiation | & E-Magnetics | EE] Output Data | K% Graphs | B3 Senstivity | 


Drive: Temperatures: E Performance Tests: 
Shaft Speed: [8000 Į Armature Winding Temperature: [40 Single operating points: 
Line Curent Definition: Magnet Temperature: [20 Open Circ 
O Peak 
@RMS 
ORMS Curent Density 2 
Feo] J jen Circuit: 
Peak Curent: [70.71 ] Shaft Temperature: |20 
za lv] Back EMF 
] erature: |20 
RMS Current: |50 | = FA Cogging Torque 
RMS Current Density: [2,847 z Electromagnetic Forces 
DC Bus Voltage: |72 | wis 
m~ On Load: 
Phase Advance [elec deg]: |0 | E-Magnetic + Thermal Coupling: Torque 
Dive E-Magnetic + Thermal Coupling: i 
Drive Type: @ No coupling (default) Torque Speed Curve 
@ Defined Currents (Defaut) O E-Magnetic Losses + Thermal Demagnetization 
O Calculated Currents O EMagnetic — Thermal Temperatures Electromagnetic Forces 
O terate to Converged Solution 
Os Mode Parameters: 
ine 
ee Skew: Sef and Mutual inductances 
Square Skew Type 
O Custom @ None (defaut) ° Transient 
Passive Generator 8 a 1 Sudden short-circuit 
Winding Connection: Rotor Stresses: 
© 
@ Star Connection (default) Centrifugal Forces 
© Detta Connection 
Magnetisation zl neat 
O Parallel 
@ Radial aes 
CO Halbach Continuous Ring Aray 


CO Halbach Sinusoidal Aray 


Figure 12. The Calculation menu where all the simulation conditions must be established. 


In this panel, the drive-related data must be set, such as the imposed speed of the 
machine, the RMS value of the supply currents, and the DC bus voltage. Additionally, 
here must be selected the type of drive, the connection of the stator windings, and the 
magnetization direction of the permanent magnets. 

By using the selection possibilities in this panel, the users can choose which perfor- 
mance tests besides the ones mandatorily computed have to be performed by the program. 
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ANSYS Motor-CAD can deliver the variations of the back EMF, cogging torque, and elec- 
tromagnetic forces (both in open-circuit and on-load conditions). There is the possibility 
to obtain the torque versus speed (mechanical) characteristics of the analyzed electrical 
machines, too. 

The simulations can be started from this panel by pressing the Solve E-Magnetic 
Model button. 


4. Simulation Results 


All the simulations to be next presented were performed at the rated conditions of the 
surface PMSM: 


e §=8000 r/min; 
50 A RMS currents; 
e 72V DC bus voltage. 


The ANSYS Motor-CAD software can perform simulations both on the entire machine 
structure or only on a part of it, considering the symmetry conditions. Despite that the FEA 
is perfectly performed at the default settings of the options concerning the solver, there 
are numerous possibilities of settings targeting the mesh fineness, the number of points in 
which the quantities of interest are computed, etc. 

Upon the performed FEA, a great variety of results were obtained. Here only those of 
major interest for the study to be performed will be detailed. 

The basic results are gained upon the accomplished FEA of the modeled electrical 
machine and can be visualized in the E-Magnetics panel of the program. For example, the 
considered regions together with the automatically generated mesh can be pictured, as can 
be seen in Figure 13. 


Figure 13. The different regions of the analyzed PMSM and the mesh generated by the ANSYS 
Motor-CAD program. 


In the same panel, the color map of four different quantities (flux density, vector 
potential, current density, and eddy current density) together with the magnetic flux lines 
can be visualized. In Figure 14, the color maps of the flux density and the flux lines for 
three distinct positions during a complete 360-electrical-degree revolution of the rotor 
are provided. 
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(b) (c) 


Figure 14. The color map of the magnetic flux density and the plotted flux lines at three different 


mechanical angle displacements of the rotor. (a) 0°; (b) 10°; (c) 20°. 


The ANSYS Motor-CAD software can generate an animation by putting together slide- 


by-slide the results obtained at different rotor positions and can save this in an animated 
gif-type file. 


A second possibility to analyze the obtained results is the evaluation of a great number 


of numeric results by accessing the Output Data panel. Here the following results categories 
(arranged in panes) can be evidenced: 


Drive: data concerning the back EMF, supply voltages and currents, d-q axes induc- 
tances and currents, machine constants, short circuit parameters, etc.; 

E-Magnetics: the main parameters of the analyzed electrical machines are provided 
here (see Figure 15), such as the electromagnetic power, torque, cogging torque, power 
factor, system efficiency, load angle, cogging period, base frequency, etc.; 

Phasor Diagram; 

Equivalent Circuit: shows all the parameters of the most complete equivalent circuit; 
Flux densities in different critical parts of the machine; 

Losses in all the parts of the machine; 

Winding: provides the parameters and current densities of all the conductors, phase 
resistance, slot filling factors, copper volume, etc.; 

Miscellaneous: data about the FEA mesh generated by the program; 

Materials: the properties and masses of all the used materials. 


As can be seen in Figure 15 and Table 1, the designed surface-mounted PMSM achieved 


all its imposed rated values (see Section 2.2). 


Table 1. The imposed rated values of the PMSM and the obtained ones via simulations. 


Item Rated Value Resulted from 


Simulation 
DC bus voltage 72V 72 V (imposed) 
Rated power 4kW 4.0596 kW 
Rated speed 8000 r/min 8000 r/min (imposed) 
Rated average torque 4.77 Nm 4.82 N-m 


The obtained power factor and efficiency of the PMSM are 0.975 and 95.36%, respec- 


tively (see Figure 15). Both are excellent values for such a machine having the given power. 
The 70.71 A phase current is under the 120 A maximum value of the power converter 
supplying the PMSM. 
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© Geometry | FJ winding | YA input Data | Hf} Calculation | £ E-Magnetics |=2] Output Data | Graphs | EG Sensitivity | 
Ay Drive [Ø EMagnetics | 2 Phasor Diagram | $8 Equivalent Crout | E Fux Densities | 3k Losses | [F] Winding | A Miscellaneous | 4 Materials | 
Variable Value Units a Variable Value Units a 
Maximum torque possible (DQ) 4,6967 Nm Flux Linkage D (Q axis current) 14,7445 mVs 
(For Phase Advance of 2,646 EDeg) 
Average torque (virtual work) 48516 Nm Flux Linkage Q (Q axis current) 3.39117 mVs 
Average torque (loop torque) 4.8292 Nm Flux linkage D (On load) 15,2883 mVs 
Torque Ripple (MsVw) 2.1162 Nm Flux linkage Q (On load) 3.45286 mVs 
Torque Ripple (MsVw) [%] 43,591 L3 — 
Cogging Torque Ripple (Ce) 1.3313 Nm Torque Constant (Kt) 0,0686552 Nm/A 
Cogging Torque Ripple (Vw) 1,5984 Nm Motor Constant (Km) 0,870424 Nm/(Watts^0.5) 
Speed limit for constant torque 8684.3 mpm Back EMF Constant (Ke) 0,0866106 Vs/Rad 
(For Phase Advance of 0 EDeg) 
No load speed 7938.4 mpm Back EMF Constant (Ke) fundamental) 0,0795389 Vs/Rad 
Speed limit for zero q axis current 1£009 pm — 
— Stall Curent 8679.79 Amps 
Hlectromagnetic Power 4067 Watts Stall Torque 595,913 Nm 
Input Power 4101.1 Watts = 
Total Losses (on load) 190,54 Watts Cogging Period 20 MDeg 
Output Power 3910.5 Watts Cogging Frequency 2400 Hz 
System Efficiency 95.354 % Fundamental Frequency 400 Hz 
= Mechanical Frequency 133,333 Hz 
Shaft Torque 4.6679 Nm Optimum Skewing Angle 20 MDeg 
Power Factor [Waveform] (lagging) 0.97579 Magnetic Symmetry Factor 1 
(User specified not to use symmetry) 
Power Factor Angle [Waveform] 12,634 EDeg Magnetic Axial Length (Slice 1) 90 mm 
Power Factor [THD] 0.9603 
Power Factor [Phasor] lagging) 0.97815 
Power Factor Angle [Phasor] 12 EDeg 
Load Angle [Phasor] 12,268 EDeg 
Phase Terminal Voltage (ms) [Phasor] 28,45 Volts 
Rotor Inertia 0,00031532 kgm? 
Shaft Inertia 4.4B69E-005 wr v EEA v 


Figure 15. The numerical simulation results provided in the E-Magnetics pane of the ANSYS Motor- 


CAD software. 


The third category of results is displayed graphically. A great variety of plots can be 
visualized from the Graphs panel of the ANSYS Motor-CAD software. The most eloquent 
ones are given in Figures 16 and 17, where the phase currents, terminal voltages, back 
EMFs, electromagnetic torque, cogging torque (computed via two different methods), and 
airgap flux density are plotted versus the rotor position. 
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Figure 16. The variation against rotor angular position of the quantities of most interest. (a) Phase 
currents; (b) terminal voltages; (c) back EMF; (d) developed torque. 
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Figure 17. The waveforms of the quantities of most interest. (a) Cogging torque; (b) airgap 
flux density. 


As the phase currents were imposed during the simulations, they have a perfect 
sinusoidal shape, and their peak value is 70.71 A, as it was forced (see this in Figure 12). 
The maximum values of the terminal voltages do not exceed the 72 V DC bus voltage. 
The back EMF produced in the generator regime is not perfectly sinusoidal, as could be 
expected. Its harmonic content can be extracted from the Harmonics window of the Graphs 
panel (see Figure 18), where the great 5th harmonic content can be easily observed. The 
generated torque has a significant ripple (equal to 43.5% upon the result given in Figure 15). 
The cogging torque is small upon the result obtained via both computation methods, based 
on the co-energy (CE) and virtual works (VWs). 
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Figure 18. The harmonic content of the back EMF waveform obtained by using ANSYS Motor-CAD 
software. 


There are some possibilities to manipulate the graphical results provided by the 
ANSYS Motor-CAD program from its Graph Editor. Various features of the plots can be 
set from here, and the final figure can be exported both in bitmap and metafile formats. 
The plotted numerical values arranged in tables can be visualized and exported in diverse 
formats for future processing in other programs. 

The authors preferred to export the results obtained through simulations to MATLAB 
via the XML file format and to take advantage of the numerous graphical processing 
possibilities existing in this very popular computing platform [68]. 
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5. Comparison of Four PMSM Types 


By using the advanced simulation features of the ANSYS Motor-CAD software, the 
comparison of four typical PMSM topologies (with surface mount, interior, spoke-type, 
and V-shaped permanent magnets) was performed. All of them have the same imposed 
rated values (given in Section 2.2), stator (concerning both their iron core and windings), 
and airgap. The rotors (shown in Figure 19) have the same permanent magnet volume. 


“NAN a L 
\ AA IN Th v 


EJ (c) A 
Figure 19. The four considered PMSM rotor topologies with different placements of the permanent 
magnets. (a) Surface mount; (b) interior; (c) spoke-type; (d) V-shaped. 


5.1. Building up the Models of the Alternative PMSMs in ANSYS Motor-CAD 


As in the ANSYS Motor-CAD software the size settings and the permanent magnet 
mass and volume computations can be performed very easily, the construction of the last 
three PMSM models was a very simple engineering task. The applied settings of the radial 
cross-section of the three PMSMs and the obtained structures are given in Figures 20 and 21. 


The settings of the axial cross-sections of the four PMSMs are identical. 
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alternative permanent magnet arrangements. (a) Interior; (b) spoke-type; (c) V-shaped. 


Stator Parameters Value Rotor Parameters Value Stator Parameters Value Rotor Parameters Yke |(nqaa panes ee io ac 
Slot Number 18 ||Pole Number 6 Slot Number Pole Number 6 Slot Number 18 || Pole Number 6 
Stator Lam Dia 118 || Notch Depth 0 Stator Lam Dia Rotor Lamination Inset 0 Stator Lam Dia 118 || Notch Depth 0 
Stator Bore 50 || Magnet Thickness 3,43 Stator Bore Magnet Thickness 3.25 Stator Bore 50 _ || Magnet Layers 1 
Slot Width (Bottom) 10,0809 || Magnet Embed Depth 05 Slot Width (Bottom) 10,0809 || Magnet Depth 143 Slot Width (Bottom) 10,0809 |/L1 Magnet Thickness 25 
Slot Width (Top) 5.9434 | Magnet Arc [ED] 113.4 | | Slot Width (Top) 5.9434. | Magnet Opening o | EE SHA ee 92 
o a heee F AR E ‘a Slot Depth 24 — |/L1 Magnet Shift 0 

| Slot Comer Radius 1 |/L1 Magnet Bar Width | 9.29 
Slot Comer Radius 1 || Airgap 0.65 Slot Comer Radius Magnet Segments 1 aa on 1 (eugene i 
Tooth Tip Depth 1 Banding Thickness 0 Tooth Tip Depth Airgap 0.65 ya 2 |L Pole V Ange m 
Slot Opening 2 ||Shaft Dia 25 See Se 7 Tooth Tip Angle 0 ||L1 Outer Extension 0.1 
Tooth Tip Angle 0 || Shaft Hole Diameter 0 Tooth Tip Angle Shaft Dia 15 Sleeve Thickness 0 — ||L1 Inner Extension 01 
eee Tai o Sleeve Thickness Shaft Hole Diameter 0 L1 Magnet Post 05 
L1 Magnet Segments 1 
L1 Mag Gap Inner 0 
L1 Mag Gap Outer 0 
L1 Layer Offset Angle 0 
Airgap 0.65 
Banding Thickness 0 
Shaft Dia 25 
Shaft Hole Diameter 0 
(a) (b) (c) 
Figure 20. The settings of the radial cross-section of the considered PMSMs with the different 
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(b) (c) 


Figure 21. The structures generated by ANSYS Motor-CAD of the alternative considered PMSMs 
having diverse permanent magnet placements. (a) Interior; (b) spoke-type; (c) V-shaped. 


5.2. The Comparison of the Obtained Results 


The simulation in ANSYS Motor-CAD of the four PMSMs was performed at identical 
settings (see them at the beginning of Section 4). 

The extraction of the multitude of results after the performed simulations enabled 
the comparison of the four PMSM topologies from different points of view. First, the total 
masses and that of each main component can be compared. In Table 2 the masses of the 
main machine components (stator and rotor iron cores, permanent magnets, windings, and 
the shaft) are given. By adding these, the total masses of the PMSMs without the housings 
and machine auxiliaries can be obtained. 


Table 2. The masses of the four PMSM components and their total masses. 


Permanent Magnet Topology 


Mass [kg] Surface Interior Spoke-Type V-Shape 
Stator core 3.736 3.736 3.736 3.736 
Rotor core 0.6537 0.7301 0.9116 0.7274 
Permanent magnets 0.1881 0.1881 0.1882 0.1881 
Windings 1.846 1.846 1.846 1.846 
Shaft 0.5743 0.5743 0.5743 0.5743 
Total mass 6.9981 7.0745 7.2561 7.0718 


Since the stators of the four PMSMs are identical, and the permanent magnets are of 
the same volume, the rotor iron core mass makes the difference. In any case, the four total 
masses are very close, and the least is that of the surface-mounted PMSM. 

In the above table and in all which follows, the best values are highlighted in green, 
while the worst are in red. 

The assessment of the main output characteristics is essential in such a comparative 
study. Therefore, eight key features of the PMSMs in the study (output power, mean torque, 
torque ripple, peak-to-peak value of the cogging torque, efficiency, power factor, power 
density, and torque density) were included in Table 3. 

As can be seen, the surface-mounted PMSM has the best characteristics in almost all 
the considered categories, while that having V-shape permanent magnets has the worst 
ones. It can be concluded that under the given conditions (identical stators, external and 
internal dimensions of the rotor iron cores, and permanent magnet masses), the topologies 
considered to be compared with the surface permanent magnet variant cannot reach the 
imposed rated speed (8,000 r/min), nor the rated output power (4 kW), nor the rated torque 
(4.77 N-m). 
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Table 3. The main characteristics of the four PMSMs. 
Permanent Magnet Topolo 
C haraci eriste Surface Interior roepe i V-Shape 
Output power [W] 4067 3516.6 3520.9 3129.4 
Mean torque [N-m] 4.851 4.192 4.193 3.731 
Torque ripple [%] 43.87 83.5 56.16 97.78 
Cogging torque peak-to-peak [N-m] 1.58 2.56 1.28 4.16 
Efficiency [%] 95.354 92.912 94.081 92.289 
Power factor [-] 0.975 0.851 0.89 0.853 
Power density [W/kg] 581.15 497.08 485.23 442.52 
Torque density [N-m/kg] 0.69 0.59 0.58 0.52 


Next, the variation against the angular displacement of some basic characteristics of 
the four considered PMSM topologies will be studied. 

In Figure 22, the variation of the generated torque in all four cases is provided. For a 
better comparison of the PMSMs, the main torque features (the minimum, maximum, and 
mean values, and the torque ripple) are given in Table 4. 
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Figure 22. The torque versus angular displacement plots for the four PMSMs. 


Table 4. The main torque features of the considered PMSMs. 


Torque Features Permanent Magnet Topology 


Surface Interior Spoke-Type V-Shape 
Tmin [N-m] 3.55 1.97 2.76 1.77 
Tmax [N-m] 5.68 5.47 512 5.41 
Tmean [N-m] 4.85 4.19 4.19 3.73 
Trippie [%] 43.87 83.5 56.16 97.78 


—___ipple "sR EEE 


As is evident, the surface-mounted PMSM has the best torque development capability, 
due to the very close placement of the permanent magnets to the airgap of the machine. 
Meanwhile, the PMSM having V-shaped permanent magnets is the worst topology, also 
from this point of view. 

In the case of all permanent magnet machines, the cogging torque is also a signifi- 
cant characteristic. Therefore, the variations of this physical quantity versus the angular 
displacement are shown in Figure 23. 
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Figure 23. The cogging torque plotted against the angular displacement for the four PMSMs. 


The peak-to-peak values of the cogging torque of the four PMSMs are significant 
comparative issues (consider Table 3). As can be seen, the lowest cogging torque is devel- 
oped by the spoke-type PMSM because in this topology the permanent magnets are deeply 
buried into the rotor iron core far from the rotor external surface and airgap. In the case of 
the V-type PMSM, the permanent magnets strongly concentrate the magnetic flux into the 
airgap, and consequentially also cause great cogging torques. 

As is evidenced, again the surface-mounted PMSM has very good behavior based on 
this consideration, as well. 

PMSMs are often used also in generator regime; therefore, the shapes of the produced 
EMF (given in Figure 24) are important issues in a such comparative study. 
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Figure 24. The back EMF generated in the case of the four PMSMs. 


As could be expected and can be observed in Figure 24, the generated waveforms are 
far from the ideal sinusoidal shapes. Hence, their harmonic content is of major concern. 
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These characteristics are again easy to extract from the results provided by the ANSYS 
Motor-CAD software. The results of the harmonic analysis were plotted by using the 
graphical potential of the MATLAB program as in the previous cases. The obtained results 
are given in Figure 25. In such cases, the total harmonic distortion (THD) of the waveforms 
is also of important relevance; therefore, their values were written on each plot from 
the figure. 
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Figure 25. The harmonic contents of the back EMF waveforms given in Figure 24. 


The PMSM having spoke-type permanent magnets has the lowest harmonic content, 
but its THD value (0.2) is only a bit smaller than the THD obtained for the surface PMSM. 
These results were expectable, since it was proved in numerous papers that these two types 
of PMSM can generate the closest EMF waveform to the perfect sinusoidal ones [69]. 


6. Discussion 


PMSMs are essential components of modern vehicles, being used both for traction 
purposes in electrical and hybrid cars and for actuating several car auxiliaries. Therefore, 
their design is an important issue and requires advanced software tools. Such a program, 
used by the authors in designing and analyzing PMSMs;, is the ANSYS Motor-CAD software. 

The simulations performed by using it enabled the authors to check the pre-sizing of 
the designed PMSM and to verify if all the imposed designed data were fulfilled. By apply- 
ing the multiple PMSM structure templates offered by the ANSYS Motor-CAD program, 
it was possible to easily perform an in-detail comparison of four rotor structures for the 
identical stator of the PMSM and to define the best structure among them. 

ANSYS Motor-CAD can also be useful in the future for improving the structure of the 
proposed PMSM by assisting in an effective optimization process. 


7. Conclusions 


By using the ANSYS Motor-CAD program in our concrete design work, we concluded 
that it was an effective tool in our hands. By using the built-in templates, the models of 
all the PMSMs taken into the study were easily achieved. Hence, the consistency of the 
pre-computed sizes of the machines was easily checked, and the CAD drawings of the 
main components could be exported for further studies or manufacturing. 

The program has a huge database of materials, where the designers probably can find 
exactly the materials proposed during the pre-sizing of the machine. If not, they can create 
their own material databases even for very specific materials. 
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The designed winding diagram could be implemented in ANSYS Motor-CAD straight- 
forwardly by setting only a few parameters. The program can display diverse diagrams 
of the modeled winding, also including the ideal displacement of the wires in the slots. 
This is one of the most spectacular features of the program since these computations are 
complicated due to the irregular shape of the slots. By the graphical illustration of the slot 
with the inserted wires and insulators, the users can easily decide on the correctness of the 
design and can have information about the actual filling factor of the winding. 

The drive systems of the simulated electrical machines can be easily configured by 
some simple selections, avoiding the use of complicated-to-build external circuits with 
adequate commands. 

By using both analytical computations and FEA, ANSYS Motor-CAD can deliver a 
great variety of results (both in graphical forms and as numerical values): the magnetic flux 
distribution in the machine, the flux densities in diverse areas, the masses of all distinct 
components, etc. Numerous output parameters (such as developed torque, cogging torque, 
power, efficiency, power factor, etc.) can be easily visualized and thus checked by the 
designing engineers. 

If needed, the correction of the design can be made rapidly until all the required 
machine parameters are achieved. 

During working with this software product, we also found some limitations of it. 
As it enables only a template-based model construction, electrical machines with special 
topologies (for example, those having modular construction or transverse flux machines), 
not included in the offered templates, cannot be studied. Despite the valuable simulation 
results obtained by using the default FEA settings, such a program should enable more 
possibilities for setting the mesh and the used solver. 

In overall conclusion concerning the used software tool: ANSYS Motor-CAD was very 
worthwhile in our design work. By applying it, we could rapidly check our designs and 
perform needed changes to improve them. 

It can be also stated that all the simulated results showed that the surface-mounted 
PMSM intended to drive the scroll-type compressor of a car HVAC system meets all 
the set design expectations. Moreover, the comparison study’s findings all underscored 
the advantages of this PMSM variant over the other three topologies that were taken 
into consideration. 
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